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Avoided band crossings have been studied theoretically and it has been shown that they can 
provide a tunning of the metal-insulator transition. Here we present an experimental example of 
an avoided band crossing for a classical undulatory system: torsional waves in locally periodic rods. 
To excite and detect the torsional waves, an electromagnetic-acoustic transducer for low-frequencies 
that we have recently developed, is used. Calculations performed using the transfer matrix method 
agree with the experimental measurements. In the observed avoided band crossing one level, which 
is a border-induced bulk level, moves from one band to the next. 

PACS numbers: 43.20.+g, 71.20.-b, 46.40.-f, 62.30.+d 



In a seminal paper, Wigner and von Neumann lj stud- 
ied the repulsion of the discrete energy levels of quan- 
tum systems when an external parameter e is varied. For 
states of the same symmetry they obtained what is called 
avoided level crossing 0. In general, however, the spec- 
trum of a quantum system is not of the discrete type. For 
instance, periodic systems show a band spectrum with 
extended wave amplitudes according to Bloch's theorem. 
If the coupling between these bands is neglected we ar- 
rive to one-band models, such as the Anderson [3| or 
the Harper models 0, |]| . When an external parameter is 
varied in these systems, some of the bands move up in en- 
ergy and others move down, giving rise to band crossings. 
At these values of the parameter, however, the coupling 
among bands can no longer be neglected. This coupling 
gives rise to avoided band crossing (ABC), which is a gen- 
eralization of the avoided level crossing mentioned above. 

The ABC has been recently studied theoretically in 
the context of the kicked Harper model Q and some un- 
expected results were found. The ABCs can provide a 
tuning of the metal-insulator transition and changes in 
the localization lengths, among other properties. 

Although the ABCs have been studied in quantum- 
mechanical systems, they can also be obtained for classi- 
cal undulatory systems such as electromagnetic, elastic or 
acoustic systems. Since we have recently dealt with elas- 
tic waves in rods with locally periodic structures, where 
wave amplitudes were measured and a band structure 
emerges |7|, we are now able to study ABC in elastic 
vibrations. It is the purpose of this letter to show ex- 
perimentally for the first time the existence of an ABC. 
Besides their intrinsic interest, the results presented here 
are relevant for vibration isolation || and in the design 
of very narrow passband filters They could also be 
of importance in applications of acousto-optic fiber de- 
vices |l0| and in the design of many aerospace and marine 
vehicles 

The rods of finite length L analyzed here consist of 
N identical unit cells of length I, as shown in Fig. ^ 
Each cell is formed by three cylinders, one of length / — e, 
radius R and cross section area S, and two cylinders of 
length e/2, radius r and cross section area s. In this way 





FIG. 1: Geometry of the rod. The length of the unit cell is I 
and e is the width of the notch. Here S and s are the cross 
section areas of the rod and notches, respectively. 



a rod with N obstacles is formed and e is the external 
parameter to be varied. We assume R « L and R « 
A = 2n/k, where A is the normal- mode wavelength, so 
the system behaves indeed as one dimensional. Torsional 
waves in this rod then obey the wave equation with a 
phase velocity y/G / p, where G is the shear modulus and 
p is the density. The normal-mode spectrum of the rod is 
pure point but, as N increases, a band structure emerges 
[jj] . 

A diagram of the experimental setup is given in Fig. 
The signal of an oscillator (Stanford Research Systems, 




FIG. 2: Block diagram of the experimental setup. The 
EMATs (exciter, detector and monitor) are shown shaded. 
The dotted line indicates mechanical contact between the cur- 
sor (filled arrow) and the scanning detector. The schematic 
diagram of the EM AT is shown in the left lower corner. 



Function Generator DS345) is sent to a power ampli- 
fier (Krhon-Hite 7500) and then to the electromagnetic- 
acoustic transducer (EMAT) exciter. This EMAT, as 
well as the one used as detector, consists of a coil and 
a magnet as shown in the left lower corner of Fig. [2] 
They are configured to excite or detect torsional waves 
without any mechanical contact with the rod. The sig- 
nal from the EMAT scanning detector is amplified by 
a high-impedance amplifier. To monitor the resonances 
the signal is sent to the lock-in amplifier (EG&G PARC 
128A). The latter converts the AC signal to a DC voltage, 
which is digitized by a CAM AC, whose output is sent to 
a PC. The DC voltage is proportional to the wave am- 
plitude when exciting a normal mode of the rod. The 
reference signal for the lock-in was taken from the oscil- 
lator. The position z of the scanning detector along the 
rod axis is measured by mechanically coupling it to a cur- 
sor in contact with a nichrome wire, depicted by the bold 
line of Fig. [21 The signal of the voltage divider is then 
sent to the CAMAC and finally to the PC. The scanning 
detector as well as the cursor are moved with a motor 
(not shown) controlled by the PC and the CAMAC. The 
normal-mode wave amplitudes are then scanned along 
the rod. We should remark, however, that the resonant 
frequencies depend on temperature. Instead of control- 
ling the latter, we changed the frequency of the exciter to 
follow the resonance. To do this, a third EMAT (denoted 
MONITOR in Fig. GJl is used to maintain the phase of 
the wave constant. The change in frequency is typically 
of the order of 0.1%, which is therefore the experimental 
error in the frequency measurements. 

The theoretical normal-mode frequencies and their cor- 
responding wave amplitudes were computed using the 
transfer matrix method 7]. In order to do this, we ex- 
press the wave amplitude ipi in cylinder i as 



i/>i(z) = A i e ik ( z - Zi ~ 1 '> + Bie-^*-**-^, 



(1) 



where Zi-i < z < Zi, i — 1, 2, . . . , 2N + 1, and impose the 
following approximate boundary conditions at the points 



where the radius changes: 



dz 
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Here Si equals either S or s according to the value of 
i. We shall consider the case of free ends, so Neumann 
boundary conditions hold at zq = or Z2N+1 = L. More 
details on the transfer matrix method for torsional waves 
in rods can be found in Ref. 0. 

The spectra of the rods depend on e/l, e being the ex- 
ternal parameter. When e/l ^ 0, 1 a band structure ap- 
pears and the bands move up or down as e/l varies. The 
coupling of the bands then leads to avoided band cross- 
ings. A typical ABC is shown in Fig. [3] for the second 
and third bands. The experimental results for e/l = 0.4, 
0.5 and 0.6 correspond to the diamonds in the same fig- 
ure. A good agreement is obtained between theory and 




FIG. 3: Resonant frequencies as a function of e/l for a 
locally periodic aluminum bar with L = 120.0 cm, R = 
6.45 ± 0.05 mm, r = 3.2 ± 0.05 mm, and N = 12. The 
experimental values (diamonds) fit rather well with the nu- 
merical ones (continuous lines). The c-level corresponds to 
n = m — 24 nodes. 



experiment. This is remarkable since no fitting param- 
eter was used. Also, from the same figure it is possible 
to see that one level, which we shall call c, moves from 
one band to the other as the ABC takes place. We will 
discuss the nature of this level below. 

In Fig. 0] we present some experimental wave ampli- 
tudes with number of nodes n — 21, . . . , 27 for e/l = 0.4 
and to = 21, . . . , 27 for e/l = 0.6; they belong to the sec- 
ond and third bands, respectively. These are the wave 
amplitudes at the left- and right-hand side of the ABC, 
which occurs at e/l = 0.5, and correspond to the nor- 
mal mode frequencies (diamonds) of Fig. |21 The results 
calculated with the transfer matrix method arc given in 
Fig.|SJ A nice agreement with the experimental wave am- 
plitudes is obtained once we take into account the fact 
that at the position of the notches the detector lies far- 
ther away than when it is located above the cylinders of 
radius R. In our case, the experimental signal is reduced 
by a factor of the order of 8, so the theoretical values at 
the notches were divided by this factor. 

As observed in these figures each wave amplitude in 
the second band is interchanged with one belonging to 
the third band located symmetrically in frequency with 
respect to the c-level, which corresponds to n — 24 in 
Figs. 0] and [3J For example, the wave amplitude for 
e/l = 0.4 and n = 22 is very similar to that with to = 26 
for e/l = 0.6. They are, however, different because they 
do not have the same number of nodes. But the differ- 
ences are negligible since they occur at points where the 
envelope of the wave amplitudes is negligible and a node 
bifurcates into three nodes. 

To quantify how much the experimental wave ampli- 
tudes of the left-hand side of Fig.0]mix with those at the 
right-hand side of the same figure, we plot in Fig. the 
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FIG. 4: Experimental torsional wave amplitudes (in arbitrary 
units) for an aluminum rod with e/l = 0.4 (left column) and 
e/l = 0.6 (right column). The amplitudes are not measured 
at the right-hand end of the rod since the signal obtained with 
the scanning detector is distorted by the exciter. 
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FIG. 5: Torsional wave amplitudes belonging to the second 
and third bands obtained with the transfer matrix method. 
The normal modes for e/l = 0.4,0.6 (left and right columns) 
are labeled by n = 21, . . . , 27 and m — 21, . . . , 27, respec- 
tively. 



absolute value of the overlap 
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between the amplitudes ip n (e/l = 0.4) and ip m (e/l = 
0.6) with n,m — 21,..., 27. The overlap is larger for 
wave functions that are symmetrical with respect to the 
c-level while it is negligible for other levels within the two 
bands involved in the ABC. This yields an anti-diagonal 
overlap matrix. The theoretical overlaps corresponding 
to the same wave amplitudes are given in Fig.HJ). Thus, 
experiment and theory agree very well and the overlap 
between theoretical and experimental wavefunctions is 
always larger than 0.84. These anti-diagonal correlations 
only appear when an ABC is present. To show this, in 
Fig. EJ; we present a numerical calculation of O nm for 
wave amplitudes at e/l — 0.4 and e/l = 0.45, since in 
this range of e/l no ABC is present. The overlap matrix 
is now concentrated along the diagonal since the wave 
amplitudes remain in its position within the band. 
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FIG. 6: Overlap O nm between wave amplitudes at e/l = 0.4 
and e/l = 0.6; a) experimental results, b) numerical results. 
In c) is shown the overlap between theoretical wave ampli- 
tudes at e/l = 0.35 and e/l = 0.4. 



So what is the nature of the c-level? It cannot be 
a Bloch-type state since at the ABC this level moves 
from one band to a neighboring one. Its existence must 
therefore be due to the fact that the system is finite so 
level c is really a border state. To corroborate this we 
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FIG. 7: Localization length of wavefunction n = 24 as a func- 
tion of x f° r e /l = 0.4. The dashed line indicates the length 
L of the rod. The insets at the left and right upper corners 
show the wave amplitudes for x = —0.03 and x — 0.03, re- 
spectively. 
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made the following transformation: z% — > Z{ + x e /2, with 
i = 1, . . . , 2N, which alters the two cylinders at the ex- 
tremes, the system remaining otherwise locally periodic. 
The c-level is then localized. The localization length as a 
function of x is given in Fig. Q We see that when \ — 0; 
the case corresponding to Figs. 0] and [5J the localization 
length is much larger than L, so the c-level appears to be 
an extended state, such as the one labeled n = m = 24. 
The c-level is actually a border- induced bulk state . 

In conclusion, we have measured for the first time an 
avoided band crossing, in this case for torsional vibrations 
of elastic rods. Since we deal with a finite rod a level 
shifts from a given band to the next one. This state is 
an example of a border-induced bulk mode, as predicted 
in Ref. [12j. We have obtained a similar behavior for 
compressional waves and it is also expected for light and 
microwaves [Kj-Hij. 
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